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ABSTRACT 

We present a model using the evolution of the stellar population in a starburst galaxy 
to predict the crystallinity of the silicates in the interstellar medium of this galaxy. We 
take into account dust production in stellar ejecta, and amorphisation and destruction 
in the interstellar medium and find that a detectable amount of crystalline silicates 
may be formed, particularly at high star formation rates, and in case supernovae are 
efficient dust producers. We discuss the effect of dust destruction and amorphisation by 
supernovae, and the effect of a low dust-production efficiency by supernovae, and find 
that when taking this into account, crystallinity in the interstellar medium becomes 
hard to detect. Levels of 6.5-13% crystallinity in the interstellar medium of starburst 
galaxies have been observed and thus we conclude that not all these crystalline silicates 
can be of stellar origin, and an additional source of crystalline silicates associated with 
the Active Galactic Nucleus must be present. 
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1 INTRODUCTION 

Silicates are among the most commonly-found dust species 
in the interstellar medium (ISM) of galaxies. Their presence 
is established through the detection of the mid-infrared res- 
onances due to the Si-O stretching and the 0-Si-O bending 
mode at 9.7 and 18 jim res pectively. In galaxie s , these bands 



star formation. The silicates observed around Asymptotic 
Giant Branch (AGB) stars can have significant crystalline 
fractions, in parti cular for the high mass-loss rate OH/IR 
stars (up to ~20%ISvlvester et al.ll 19991 : iKemper et "alll200ll : 



are s een in absorption (iGillett et al 



19751: ISturm et all 



2000l). as well as in emission (jSiebenmorgen et al.l l2005l : 



Hao et al.|[2005D . Most of these silicates show the broad 
resonances characteristic of amorphous silicates, i.e., sili- 
cates showing a large degree of lattice defects, and it is 
generally assumed that silicates in the ISM are predomi- 
nantly amorphous. In particular, the degree of crystallinity 
X, defined as the mass fraction of silicates that is crys- 
talline, X = Mx/{Mx + Ma), in the Galactic diffus e ISM 
is fo und to be of the order of 1% (|Kemper et al.l l200ll . 
I2OO5I ). This contrasts sharply with the much higher de- 
gree of crystallinity seen in silicates in the circumstellar en- 
vironments of pre- and post-main-sequence stars (see e.g. 
iMolster fc Kempeill2005l . and references herein). Generally 
speaking, the galactic cycle of dust starts with its formation 
in evolved stars, followed by processing in the ISM and even- 
tually ends with incorporation in stars and planets during 
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Ide Vries et"ai]|2010l ). For lower mass-loss rate AGB stars 
such as Miras, the crystallinity is not well established, but 
is co nsistent with a value that does not vary with mass- loss 
rate (|Kemper et al.ll200ll V For more massive stars, such as 
Red Supergiants (RSGs), the crystallinity of the silicates in 
the stellar ejecta is not well known, alt hough isolated studie s 
report high crystalline fractions (e.g. iMolster et al.l 0*9991 ) . 
while low crystalline fr actions seem to be more common 
(|Verhoelst et all |2009| ). ISpoon et all (|2006l ) adopt a crys- 
tallinity of 15% for RSGs and Luminous Blue Variables. 

In contrast to the low crystallinity in the Galac- 
tic ISM, significant amounts of crystalline silicates have 
been detected in the infra red spectra of Ultr aluminous In- 
frared Galaxies (ULIRGs: ISpoon et al.|[2006l ). In a sample 
of 77 ULIRGs, 12 were found to show crystallinity, with 
crystalline-to-amorphous silicate mass ratios ranging from 
0.07 to 0.15, corresponding to crystallinities of 6.5% to 13%. 

The degree of crystallinity of a population of silicate 
grains provides a record of the processi ng history of those 
grains (see e.g. IMolster fc Kempeijl2005l V A high crystalline 
fraction points to a relatively high formation or processing 
temperature (~ 1000 — 1500 K), while a large amorphous 
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fraction indicates that the population of grai ns has under- 
gone the damaging effects of cosmic ray hits (jBringa et al.l 
|2007^. grain-grain collisions or atomic impacts in shock s 
l|Demvk et aij|200ll : iBrucato et aLll2003l : IJager et al.ll2003h : 
or that it is formed at lower temperatures. Crystalline sil- 
icates may thus form in dense circumstellar environments, 
where the vicinity of the central star provides the required 
heating. The mass-loss processes of evolved stars subse- 
quently spread these crystalline silicates into the ISM. The 
fact that the silicates in the ISM of the Milk y Way are al- 
most entirely amorphous (|Kemper et al.|[2004l 'l suggests that 
the amorphisation processes in the ISM are more impor- 
tant than the injection of fresh crystalline silicates into the 
interstellar reservoir by evolved stars. For the Galaxy, an 
amorphisatio n time scale of 40 Myr has been derived from 
observations (|Kemper et al.l |2004| . l2005l). which is cl ose to 
the experimental value of 70 Myr (|Bringa et ahllioOTl ). 

ISpoon et all (|2006l ) argue that in starburst galaxies ul- 
traviolet (UV) photons from either the Active Galactic Nu- 
clei (AGN) or massive stars are the only potential sources 
of UV photons to anneal or form the crystalline silicate 
mass observed. Crystallisation due to UV photons orig inat- 
ing from the AGN is dismissed bv lSpoon et al.l (|2006l ). be- 
cause of the observed l ack of crystalline silicates in the in- 
ner 2 pc of NGC 1068 (|jaffe et al.ll2004l '). and the fact that 
the crystalline silicates are only seen in absorpt ion, and can 
therefore not be very warm. ISpoon et al.l (|2006i l hypothesise 
that the crystalline silicates must be produced by massive 
stars originating from the starburst. In this paper, we will in- 
vestigate the viability of the build-up of crystalline silicates 
due to the starbu rst activity in ULIRGS as proposed by 
ISpoon et al.l (|2006h . and compare the crystalline fraction to 
the levels observed by these authors. The alternative forma- 
tion of crystalline silicates due to AGN activity will be the 
subject of a separate future study. This future work will be 
based on the predicted format ion of dust in the quasar wind 
rising from the accretion disk (|Elvis et al.ll200 , where con- 
ditions may be similar to those present in AGB star winds, 
a class of very efficient dust producers. Crystalline silicates 
may form in quasar winds if the conditions are right. Indeed, 
we have already observ ed crystalline silicates in the q uasar 
wind of PC 2112-^059 (|Markwick-Kemper et ai]l200'iD . 



2 THE MODEL 

2.1 Injection of mass in the ISM of starburst 
galaxies 

When a starburst occurs, a large fraction of the gas available 
in molecular clouds will be triggered into forming s tars. 

Essentially following the method described by iTinslevI 

l|l980t) , we assume that the burst of star formation occurs 
at a constant rate ip{t) — ijjo for to < t < ti, and that the 
star formation rate before to and after ti is negligible. The 
distribution of stars over the mass spectrum follows the ini- 
tial mass function (I MF) (/)(m ) = iVo m~^^+°'. We adopt the 
IMF as described bv lKroupal (|200lD . who derive a value of 
a = 1.3 for M > 1 Mq, which is not too dissimilar from the 
original Salpe ter IMF (a = 1.35 for the entire mass range; 
ISalpeteij[l955l ). About 51% of the mass in a star formation 
event is conta ined in the stars more massive than the Sun 
jKroupalbOOll '). from which we can calculate A^o by solving 



/ (t>{m)dm = No m"<^+"'dm = 0.51. (1) 

im=l Jl 

We find that A'o = 0.663 satisfies this equation. Since 
the duration of the starburst will be much shorter than the 
life time of the Sun (r©), we do not need to consider mass- 
loss for low mass stars, so the exact shape of the ISM below 
M < 1 Mq is irrelevant. Naturally, setting the upper inte- 
gration boundary to m = oo results in unrealistically high 
stellar masses, but since the IMF is a very steep power- 
law, the numerical contribution from hypothetical stars with 
M > 120 Mq is negligible. 

The time that a star spends on the Main Sequence is 
approximated by tms oc while the mass- luminosity re- 
lation states that L oc Af, with rj ~ 3.35 for stars wit h 
IMq < M < 40 Mq (e.g. iKippenhahn fc Weigei^ liggd l. 
Thus we find for the lifetime of a star on the Main Sequence 
- using solar units, e.g. m{MQ) and t{TQ) - that Tm ~ m}~^ . 

It is possible to set up an equation for the ejection rate 
E{t) at which gas is being expelled from a population o f 
stars formed during a starburst in a galaxy (|Tinslevlll980l V 
Assuming the mass-loss occurs instantaneously after a life- 
time Tm on the Main Sequence, a star formed at time {t — T,n) 
ejects its shell at time t. In order to obtain the ejection rate, 
one needs to integrate over all masses m larger than nit, 
the stellar mass corresponding to lifetime Tm = t- Thus, the 
ejection rate is 

/"OC 

E{t) = / (m - Wm)i'{t - T,n)(j>{m)dm (2) 

where Wm is the mass of the remnant after the mass loss. By 

, _ 1 
substituting Tm = rn ^ and mt — t i-i and using the initial 

mass function and our assumption that the star formation 

rate is constant at a value of ijja during the burst, the ejection 

rate E(t) can be written as 

E{t)= _^{m - Wm)ipoNom-^^+''^ dm . (3) 

J m=t l-l 

2.2 The remnant mass function 

We a dopt the strategy utilised by IDabringhausen et al.l 
(|2009l) to describe the remnant mass function, Wm- It slices 
the stars with m ^ 1 into three categories: those that 
leave white dwarfs (m < 8), those that leave neutron stars 
(8 ^ m < 25) and tho se that leave black ho l es (m ^ 25) as 
their stellar remnants. IDabringhausen et al.1 (|2009l l arrive at 
the following expressions for the remnant mass of the form 
Wm ~ fmm + Wmg , with /„, aud Wmo determined in each 
category to be: 

{0.109 m-f 0.394 for l<m<8 
1.35 for 8 < ?n< 25 (4) 

0.1m for 25 < m 

2.3 A starburst of finite duration 

In order to calculate the ejection rate, we can distinguish 
two different eras: during the starburst (0 < t < ct) and af- 
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ter the starburst [t > a). Within each era, we can subdivide 
into three diflerent regimes, correspo nding to the lifetimes 
of the stellar mass ranges described bv lPabringhausen et alj 
l|2009h . defined by the boundaries rg = 8^"''r0 = 7.5 • 
10"^ T0 and r25 = 25^~''r0 = 5.2 ■ 10"" tq. We assume 
that o" > rg. 



2.3.1 During the starburst: <t < a 
• Phase 1: t < T25: 

Only the most massive stars (m > 25) have started losing 
mass, and for these objects the remnant mass is given by 
Wm = 0.1m (Eq. lU, and thus Eq. |3] becomes 

/" 00 

E{t) = / m(l - /„)V^oiVom-<^+''>dm (5) 

J m — m± 

which simplifies to 



E{t) = -i/;oiVo(f - U) 



1 

mt 



1 - a 



(6) 



with iVo = 0.663 and = 0.1, and rut = {t/T^f^''^. 
• Phase II: r25 < t < rg: 

During this phase, two components contribute to the ejec- 
tion rate: E{t) = _Ei + Eiit), with E-i a constant value due 
to the ejection rate by the most massive stars (m > 25), 
and E2{t) an increasing component due to the ejection by 
intermediate mass stars (8 < m < 25). To obtain Ei, we 
integrate Eq. [5] between m = 25 and 00, and find that 



El = -(1 - ^)7/;oiVoY-^25'"" 



(7) 



which, for = 0.1, A^o = 0.663, and a = 1.3 results in 
El = 0.76 V-o. 

To obtain E2{t), we integrate 



E3it)= / 1 {m-Wm)ipoNom-'-^+''^dm (11) 

J m — mt — {t/ tq ) ^ ~ ^ 

with Wm ~ fwin + Wmg ■ This can be rewritten as 



E-iit) = i:oNo 



1 - a 



a 8 



-mt 



1 - .fu, 

1 - a 



a mt 



(12) 



2.3.2 After the starburst: t > a 

When a starburst of duration a has stopped, there are no 
longer any stars with m > {{t — a) /tq) 1-1 in the galaxy, 
thus constraining the upper limit of the integration to obtain 
the ejection rate. 

• Phase IV: t — a < T25 : 

Shortly after the starburst has stopped stars with m > 25 
still exist, and we can write for the ejection rate: E{t) = 
E[{t) + E2 + Esit). E2 and Eait) are given by Egs.irUland 
[T2I For Ei{t) we can write, with mt-a = {(t - cr)/r0)^/^"'', 

E[{t)= m(l - /»)i/;oAfom-''+"'dm (13) 

J m=25 

which we can rewrite as: 



E{t) = i:oNo{l - U) 



1 - a 



mt- 



-25'"" 



(14) 



• Phase V: r25 < t — (J < rg. 

In this phase, high-mass stars (m > 25) are extinct and no 
longer contribute to the ejection rate, which has now become 
E{t) = E'2(t) + Ei{t), with E3{t) again equal to the result 
of Eq. [121 For E'2 (t) we can write 



E2{t)= {m-Wm)^oNom-^^+''^dm (8) 

J m — mt 



and find 



E2{t) = VoiVo 25 \ ] - mt 1 

\ Vl — a a / Vl — a a / 

(9) 

with w,n = 1.35. 

• Phase III: rg < t < a, with a < tq: 

During this phase, the ejection rate can be written as 
E{t) = El + E2 + Esit) with. El given by Eq. [T] and E2 
obtained by substituting mt — 8 into Eq. |S] This yields 



E'2it) 



m = 8 

VoiVo 



mt- 



mt~ 



-a_ ^ Uhn 

~ a a 
a a 



(15) 



with mt-„ = ((t-cr)/r0)i/i-''. 
• Phase VI: t — a > ts 

At this late time after the starburst only low-mass stars 
with m < 8 are still in existence, and thus 



E2 = ^oA^o 25 



25 
1 - a 



Wm 

a 



8 Wm 

1 — a a 



(10) 

which, for Wm = 1.35 gives E2 = 0.47 ipo. 

The ejection rate by low mass stars (m < 8) Ez{t) can be 
obtained from Eq. [S] 



l~a 1~ fw , Wmo 
mt-a \ h 



E{t)^E'-,{t) 
1 



-mt 



1 — a a mt 

l-a / 1 ^ /ill I ™mo 1 



1 - a 



a mt 



(16) 
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Figure 1. Stellar mass ejection rate into the ISM of a starburst 
galaxy with a star formation rate -0 = 1000 Mq yr"-*^, as a func- 
tion of time after the beginning of the starburst (solid line). The 
duration of the starburst is 100 Myr. The dotted, dash-dotted 
and dashed lines indicate the contributions by high, intermediate 
and low mass stars respectively. 
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Figure 2. Dust production rate D{t) in a starburst galaxy, com- 
pared to the stellar mass ejection rate E{t), both indicated by 
solid lines. The dust-to-gas ratio d is assumed to be 0.01 in all 
ejecta. The dashed and dotted lines show the dust production 
rates in cases where supernovae produce dust at 10 or 100 times 
lower efficiencies, respectively. 



2.4 Mass ejection rates and dust production 

In Fig. [1] the resulting ejection rate E{t) throughout phases 
I- VI is shown for a starburst of 1000 Mq yr~^ and a du- 
ration of 100 Myr, as well as the ejection by high, interme- 
diate and low mass stars specifically. It is clear that during 
the burst of star formation supernovae ejecta from the most 
massive stars (m > 25) dominate the mass return to the 
ISM, but that after the starburst has ended, intermediate- 
and low-mass stars are important. 

For each of the regimes, the dust production rate D{t) is 
related to the mass ejection rate by D{t) = E{t)d, in which 
d is the dust-to-gas ratio. In the simplest scenario we assume 
that low, intermediate and high mass stars are all equally ef- 
ficient dust producers, for which we use the commonly-used 
value of d = 0.01. This gives rise to the solid line labelled 
D{t) in Fig. [2] The dust-production efficiency of supernovae 



is a hot topic in research because of its implications for the 
presence of dust at high redshift, and the consensus seems to 
be that dust production in supernova ejecta is not as efficient 
as in the outflows of lower mass stars (e.g. ISugerman et al.l 
l2006l : lFox eTallbOlOl ). We have thus also considered values 
of 0.001 and 0.0001 for d, for the ejecta of high and interme- 
diate mass stars (m > 8), all of which produce supernovae. 
The resulting D{t) is shown in Fig. [2] by the dashed and 
dotted lines respectively. 



2.5 The dust budget in the ISM of a galaxy 

In our model, we will consider two separate dust reservoirs 
to be present in a galaxy: a crystalline silicate reservoir with 
a mass Mx and an amorphous silicate reservoir with a mass 
Ma - The injection of new silicate material D{i), as described 
in Sect. 12.41 affects both reservoirs, although not in equal 
amounts. We use a;* to denote the fraction of silicate mass in 
crystalline form produced by stars. Another possible source 
of silicates is dust formation in the ISM itself, which can be 
generally expressed as FA{t) and Fy^[t) for amorphous and 
crystalline silicates respectively. 

We include the transitions between both reservoirs due 
to annealing or destruction of the lattice structure, which 
occur at a rate fci for the amorphisation of crystalline ma- 
terial and k2 for the crystallisation of amorphous silicates. 

Finally, the processes that remove material from both 
dust reservoirs include grain destruction by supernova 
shocks and cosmic ray hits, but also astration and incorpo- 
ration in planet-forming disks. We summarise these terms 
as dust destruction, occurring at rates of fca and ki for crys- 
talline and amorphous sili cates respectively. 

Following Eq. 5 from lKemper et al.l (|2004l ). these rates 
can be written as a set of coupled differential equations: 



dt 

dt ^ \^ 



fciMx + fcaMA - fesMx + Fx{t) 

■ X,)D{t) + fclMx - feMA - fc4MA + FA{t) . 

(17) 

In interstellar conditions, the crystallisation rate k2 will be 
negligible, as annealing requires the grains to be heated to 
> 1 000 K, which is unlikely to occur. We also assume that 
the destruction rates fo r both types of silic ates are equal: 
fca = fc4 = 2 X lO"'' yr"^ (|Kemper et al.ll2004l ). and that dust 
formation in the ISM can be ignored: Fy^{t) = Fxit) — 0. 
For the amorphisation rate of crystalline silica tes we take the 
rate ki — 2.5 x 10~* yr~^, in accordance with iKemper et al] 
(|2005i ). 

Eq. [T7] allows us to investigate the time dependence of 
the crystalline and amorphous silicate masses, and compare 
the r esults with the observational constraints (|Spoon et al.l 
l2006t ). Instantaneous mixing of the silicates is assumed. 



2.6 Fitting parameters 

The values for the fixed parameters and the fitting ranges 
for the free parameters required to reproduce the observed 
crystallinities are set as follows. The star formation rate ipo 
is allowed to vary between 10-1000 Mq yr~^. The total ini- 
tial dust mass in a galaxy is varied between 10^-10* Mq, 
with an initial crystalline fraction a^isM ^ 0.01. It should be 
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noted that both the initial dust mass and the star forma- 
tion rate are related to the gas mass in the galaxy, and both 
will go up for a larger gas reservoir. The crystallinity in the 
stellar ejecta a;, is taken to be 0.10-0.20. The duration of 
the starburst is allowed to run until 0.01 tq, or 100 Myr, 
which is well beyond the commonly accepted value of 5-10 
Myr for the duration of starburst, although recent work on 
dwar f galaxies suggest tha t starburst may last a few 100 
Myr l|McQuinn et al.ll2010t ). 



3 RESULTS 

We have numerically evolved the system using Eq. 1171 The 
results of a run with an initial dust mass Mx + Ma = 
10* Mq, star formation rate t/jq ~ 1000 Af© yr~^ and crys- 
tallinity I, = 0.2 for the silicates in the stellar ejecta in 
Fig. O In this case supernovae are considered equally effi- 
cient dust producers as low mass stars, with d = 0.01 over 
the entire stellar mass range. Panel a) shows the silicate 
ejection rate throughout time for a starburst of 0.01 tq in 
duration. Due to the contribution of these ejecta, the total 
silicate mass Mx -I- Ma in the galaxy increases (panel b)), 
and levels off after the end of the starburst (dashed line). 
The dash-triple-dotted line indicates the destruction time 
scale of 400 Myr, which is due to rates and k^. At this 
point the total dust mass has decreased significantly from its 
peak during the starburst. The mass Afx contained in crys- 
talline silicates (panel c)) closely follows the ejection rate, 
since the amorphisation time scale is short (40 Myr, based 
on fci , indicated with a dotte d line). The crystalli nity in the 
ISM xisM = Mx/(Mx+Ma) (|Kemper et al.ll2004l ') is plotted 
in panel d) and appears to peak before the amorphisation 
time scale has elapsed, at around 20 Myr. A significant level 
of crystallinity may build up for a brief period of time. In 
this particular example, we find levels of a;isM ~ 10 — 11%. 

The definition of crystallinity use d here is based on th e 
work presented in an earlier paper jKemper et al.ll2004h . 
and differs fr om the crysta ll ine-to -amorphous ratio A'cr/A^am 
measured bv ISpoon etal] (|2006l '). Their ratios of 0.07-0.15 
therefore translate to crystallinities xism of 6.5% - 13% 
for the interstellar silicates in these ULIRGs, experienc- 
ing starbursts. The peak crystallinity in our calculation of 
~ 10 — 11% falls within this range. 

However, a star formation rate of Tpo — 1000 M0 yr~^ 
is on the high end of the rang e seen in starburst galaxies 
(e.g. IWeedman fc Houckll2008l'). and values o f 10-100 seem 



to be more common (e.g. Kennicut3 Il998l ). When these 
lower values are used in the model calculations, the crys- 
tallinity of the silicates in the ISM builds up more slowly 
(Fig. HI), and peak s at values below the measurements by 
ISpoonet all (|2006l ). 
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Figure 3. Results from model calculations, with the initial in- 
terstellar silicate mass Mx + Ma = 10* Mq , the star formation 
rate i/jq = 1000 Mq yr"-*^, the stellar crystaUinity a:* =0.2 and 
the dust-to-gas ratio d = 0.01 for all stars. Panel a) displays the 
total injection rate of stellar dust into the ISM as a function of 
time after the beginning of the starburst. In panel b) the total sil- 
icate mass in the galaxy is shown. Panel c) shows the crystalline 
silicate mass in solar masses, and finally, in panel d) we show the 
crystallinity, defined as the crystalline silicate mass divided by 
the total silicate mass. The dashed fine indicates the end of the 
starburst ti = 10* yrs. The dashed-triple-dotted and the dotted 
lines indicate destruction and amorphisation time scales in the 
ISM respectively. 



3.1 The effect of supernovae 

As Figs. and U demonstrate, the levels of crystallinity 
reached in starburst galaxies with high (1000 Mq yr~^; 
Fig. |3ll and intermediate (> lOOM© yr"^; Fig. [l]) star for- 
mation rates are i n line with the ran ge of a; = 6.5 — 13% 
that is observed bv lSpoon etal] (|2006D . However, these cal- 
culations were performed under the assumption that the 



dust production by supernovae is efficient, with a dust-to- 
gas mass ratio d in the stellar ejecta of 0.01, independent 
of initial stell ar mass. However, when we, fo l lowin g work by 
e.g. rSugcrman et al.1 (120061 ) and IFox et al.l (|2010l ') confirm- 
ing that supernova are not efficient dust producers, consider 
a 10 or 100 times lower dust formation efficiency in the su- 
pernovae ejecta (affecting all stars with m > 8) a different 
picture arises. 



6 F. Kemper, A. J. Markwick and Paul M. Woods 



10% 


1 1 1 1 — 1 1 1 1 1 1 1 1 1 — 1 1 1 1 1 1 1 1 1 — 1 1 1 1 

/ \ ' 

/ooo \i 


>, 








B 




crys 


/ 100^ k\ 





^^^^^^^ -i^^O^^^ 





I I I I ' ' ' I I I I I I ''III I I I I III' 

10^ lO' 10** lo'' 

t(yr) 



Figure 4. Crystallinity in models with star formation rates tpo of 
1000, 100 and 10 Mq yr"! (solid lines). The other model param- 
eters are set as described in the beginning of Sect.fS] The dashed 
line indicates the end of the star burst. 

The top two panels of Fig. [S] show the amorphous and 
crystalline silicate mass in the starburst galaxy, using values 
of d = 0.01 (solid lines), 0.001 (dashed lines) and 0.0001 
(dotted lines) describing the supernova dust productivity, 
while the dust production efficiency of low mass stars (m < 
8) is kept at d = 0.01. The two sets of lines in the top panel of 
Fig. [S] arise from different initial dust masses in the ISM (see 
caption) . These plots show that the overall dust production 
dramatically decreases, and that the peak in dust mass shifts 
to a later time. The crystalline silicate mass is even 1 or 2 
orders of magnitude lower during the starburst phase, and 
is no longer dependent on the supernova dust productivity 
20 Myr after the starburst. This is due to the fact that the 
amorphisation time scale is very short, and all crystalline 
silicates present at this time are freshly produced by low 
mass stars. 

The crystallinity x, in the bottom panel of Fig.[5]shown 
for Tpo — 1000 Af0 yr~^ and an initial dust mass of 10* Mq, 
evolves differently with time for different dust production 
efficiencies in supernova, however the width of the curve 
is determined by the amorphisation rate fci. At the lowest 
efficiency we considered, d — 0.0001 the crystalline fraction 
is in fact dominated by the produce of low mass stars, which 
peaks after the end of the starburst, in this case. The peak 
crystalline fraction is also considerably lower for low dust 
production efficiencies, and for this very high star formation 
rate only a level of 4% is r eached, just below the lowest 
values of 6.5% observed bv ISpoon et all (|2006l ). After 5- 
10 Myr, the generally accepted value for the duration of a 
starburst, only small amounts of crystallinity have built up, 
even for the most favourable set of parameters. 



4 DISCUSSION 

Through constructing an ISM dust evolution model we 
have been able to explain the crystallinities observed by 
ISpoon et"aD (j2006l ) in a sample of 12 ULIRGs, out of a 
larger set of 77 analysed. The observed crystallinities range 
from 6.5 to 13%. We remark, however, that the derived crys- 
tallinity is high, due to the use of the 16.5 fj,m forsterite 
feature. This is intrinsically one of the weaker forsterite 
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Figure 5. The amorphous and crystalline dust masses in a star- 
burst over time, for three different SN dust production efficiencies: 
d = 0.01 (solid Hnes), d = 0.001 (dashed lines) and d = 0.0001 
(dotted lines). The first panel shows the amorphous silicate mass 
Ma in the ISM of a starburst galaxy, with the two curves in each 
line style representing the initial ISM dust masses of 10^ Mq 
(lower lines) and 10® Mq (upper lines). The second panel shows 
the crystalline silicate mass Mx, and the third panel the crys- 
talline fraction x = Mx / {Mx + Ma) of silicates in the ISM of a 
starburst galaxy as a function of time. For all models, Vo is set 
to 1000 Mq yr-^ 

resonances (|jager et al.lll993 ). and thus a large amount of 
forst erite is required to e xplain the optical depth in this fea- 
ture. ISpoon et al.) (|2006l ) chose this feature for their analysis 
because it was the only feature clearly observed in all 12 
objects. Using Fig. 5 and Eq. (1) from their paper, we esti- 
mate that the crystallinity could be a factor of ~3 or ~4.5 
lower, when using the 18 or 23 /im features respectively. This 
should, however, be compared to the lack of a crystalline sil- 
icate detection in the ISM of the Milky Way jKemper et al.l 
I2OO4I ). Additionally, crystallinities in the ULIRG sample can 
only be derived for all sources using the 16.5 ^m featu r e, and 
we therefore use the results presented bv lSpoon etHI (|2006l ) 
at face value. 

Using the relatively high values of 1000 Mq yr^^ for the 
star formation rate and a crystallinity of 20% in the stellar 
ejecta we achieve interstellar crystallinities of > 10%, under 
the assumption that supernovae are efficient dust produc- 
ers, with a dust-to-gas ratio d — 0.01 in their ejecta. As the 
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time-dependent plots show, detectable levels of crystallinity 
are only a temporary condition, which is consistent with 
only seeing it in part of the ULIRG sample studied, with 
the timing and the duration of detectable levels depending 
on the model parameters. In the case that supernovae are 
important dust producers, the crystallinity of the silicates 
in the ISM peaks shortly after the starburst on a time scale 
shorter than the amorphisation time scale of 40 Myr. For the 
parameters used in our fiducial model, we find that the crys- 
tallinity is highest some 10-20 Myr after the starburst and 
that crystalline silicates can make up of the silicate 

material. 

4.1 Evaluation of the input parameters 

The input parameters used are on the extreme end of the 
range observed, particularly for the star formation rate V'o, 
the crystallinity x, of the silicates in stellar ejecta and the 
(supernova) dust formation efficiency d. We also discuss the 
role of supernovae in general. 

First, for more typical star formation rates of V'o = 
10 — 100 Mq yr^^ the interstellar crystallinity is much less 
enhanced during the starburst (see Fig.|4|, making it harder 
to explain the observed values. 

Second, the interstellar crystallinity scales linearly with 
the crystallinity of stellar ejecta used, where x* = 20% 
is on the high side. More typi cal values are around 10% 
(|Kemper et all 1200 ll ). although [de Vries et all (|2010l ) have 
shown that the forsterite abundance (one of the species 
contributing to the crystallinity) is closer to 12%. Unfortu- 
nately, the crystalli nity can only be mea sured in stars with 
the densest winds (|Svlvester et al.lll999l ): crystallinities up 
to about 20% in the more numerous low d ensity winds would 
remain undetected jKemper et al.ll200ll '). Thus, 20% could 
be regarded as an upper limit to the crystallinity in stellar 
ejecta, for classes of objects where no specific measurements 
exist. 

4-. 1.1 The role of supernovae 

The role of supernovae in the mineralogical evolution of a 
galaxy remains uncertain, for a number of reasons. First, lit- 
tle is known about the mineralogy of su pernovae ejecta. The 
small number of studies available (e.g. iRho et 811120091 ) do 
not establish the crystallinity of the produced silicates. Sec- 
ond, the dust product ion efficiency d of supernovae has been 
shown to be low (e.g. ISugerman et al.ll2006l : lFox et al.ir2010l ') 
compared to lower mass stars. On the other hand, studies 
of the Gas A and Kepler supernova remnan ts suggest sig- 
nificant amounts of dus t may have formed (|Morgan et al.l 
l2003l : Id unne et al.ll2003l ). although in one of these particu- 
lar studies foreground emission may have bee n interpreted 
as co ming from the supernova remnant itself jKrause et al.l 
\2004 ). Thus, d may range from 0.0001 to 0.01. In Fig.0 the 
effect of a low dust formation efficiency in supernovae on 
the dust population in a galaxy shows that the crystallinity 
does peak at much later times and at lower values due to 
his effect. Indeed, at a supernova dust production rate of 
d = 0.0001, the dust production, and thus the crystallinity, 
is dominated by low-mass stars (m < 8, with d = 0.01), 
and by the time the crystallinity peaks these galaxies are no 
longer recognisable as starburst galaxies. 
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Figure 6. Crystalline fraction of silicates x = Mx/{Mx + Ma) 
in the ISM of a starburst galaxy as a function of time, for three 
different amorphisation rates k\ = 2.5 ■ 10~* (solid line), fci = 
2.5 ■ 10"^ (dashed line) and fci = 2.5 ■ 10"'' yr"! (dash-dotted 
line). The curves shown are for the initial dust mass of 10* Mq, 
and a star formation rate of 1000 Mq yr^-*^. 



Third, the higher supernova rate has its effect on the 
amorphisation rate. The values used for fci are derived for 
Galactic conditions but these may be significantly higher in 
starburst galaxies due to the enhanced rate of shocks and 
cosmic ray hits. We have investigated the effect of an en- 
hanced supernova rate, by increasing fci up to a factor of 
100. Increasing these rates is detrimental to the crystalline 
fraction, with an increase of a factor 100 reducing the crys- 
talline fraction in the ISM to <l%o (Fig.©. 

Finally, as ISpoon et al.l (|2006l ) argue, supernova pro- 
genitors, such as RSGs and Luminous Blue Variables, 
may produce significant amount of crys talline silicates 
(e.g. iMolster et al]|l999l : IVoors et al.lll999l '). and thus con- 
tribute to the interstellar crystallinity. However, the sub- 
sequent supernova shocks travelling through these earlier 
ejecta will destroy a large fra ction of the grains by sput- 
tering (jjones et al. I I 1994 I 19961 ). The crystalline and amor- 
phous silicate materials (assumed to be affected in equal 
amounts) will be returned to the gas phase. Thus, the ISM 
will be enriched by the atomic building blocks of silicates, 
but in gaseous form. Moreover, the high supernova rate in 
starburst galaxies will affect the dust destruction rates fcs 
(= ^4), returning even more dust to the gas-phase. 



4-1.2 Dust formation in the ISM 

The (re-) condensation of dust from the gas phase in the ISM 
will also affect the crystallinity. In our model calculations, 
we assumed that the dust formation in the ISM itself is neg- 
ligible as a sou rce of dust, hence Fa = Fx = 0. However, 
calculations bv lZhukovska et al.l (|2008l 'l show that it can be 
the most important process that contributes to the dust con- 
tent of the Milky Way. Since dust formation and dust growth 
at the low temperatures prevalent in molecular clouds yields 
only amorphous silicates, inclusion of this factor will further 
reduce the crystallinity of interstellar silicates. 
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4.2 Low crystallinities in typical starbursts 

Although our models are succe ssful in explaining the crys- 
tallinity observed in ULIRGs (|Spoon et alJl2C)0B ). we find 
that this is only the case for an extreme set of input param- 
eters. For more usual values for the star formation rate and 
the dust formation efhciency in supernovae, and when the 
effect of supernovae on grain destruction and amorphisation 
is taken into account, we find that the observed crystallinity 
of silicates cannot be explained by a starburst driven model 
only. Thus , we conclude that the high level of crystallinity 
derived by ISpoon et al.l (|2006l ) is inconsistent with the for- 
mation of crystalline silicates in stellar ejecta only. However 
their observations are also consistent with crystallinities a 
factor of ~ 3-4.5 lower than what is reported in their paper 
when different resonances are used, in which case our model 
calculations show that stellar sources could indeed explain 
the observed crystallinity, using reasonable parameters. 



5 CONCLUSION 

In order to explain the h igh crystallin e fract ion of silicates in 
ULIRGs, as reported bv lSpoon et al.l (|2006l ) we find that an 
additional source of crystalline silicates (or crystallisation of 
amorphous silicates) must be present, related to the AGN 
itself, rather than the starburst activity. In this case, the 
crystallisation will occur due to heating by (UV) photons 
from the AGN environment, rather than (massive) stars. A 
potential scenario may in clude the formati on of (crystalline) 
silicates in quasar winds (|Elvis et al.ll2002l '). In addition, fur- 
ther observational studies are useful to better establish the 
crystalline fraction of silicates in ULIRGs, and validate the 
conclusions presented here. 
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